The halo of the exotic nucleus ^^Li: a single Cooper pair 
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■ If neutrons are progressively added to a normal nucleus, the Pauli prin- 

(N 

ciple forces them into states of higher momentum. When the core becomes 

(N 
> 
m 

■ neutrons outside to form a halo, which because of its large size can have lower 

■ momentum. It is an open question how nature stabilizes such a fragile system 
and provides the glue needed to bind the halo neutrons to the core. Here we 
show that this problem is similar to that of the instability of the normal state 

^ . of an electron system at zero temperature solved by Cooper, solution which is 



X 



neutron-saturated, the nucleus expels most of the wavefunction of the last 



at the basis of BCS theory of superconductivity. By mimicking this approach 
using, aside from the bare nucleon-nucleon interaction, the long wavelength 
vibrations of the nucleus ^^Li, the paradigm of halo nuclei, as tailored glues 
of the least bound neutrons, we are able to obtain a unified and quantitative 
picture of the observed properties of ^^Li. 
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Aside from "dark matter" |]T|, atomic nuclei, little droplets made out of protons and 
neutrons 10 femtometer across [0], make up a sizable fraction of the present mass of the 
universe. Research into the structure of atomic nuclei concentrates largely on the limits of 
the nuclear stability ||^, where new physics is expected to be, in particular at the limits of 
neutron and proton number defining the so called drip lines in the chart of nuclides. The 
most exotic nuclei [Q, first produced in the laboratory only few years ago, are those that 
lie just within the drip lines, on the edges of nuclear stability. Of these, the atomic nucleus 
s^Lig, containing 3 protons and 8 neutrons, is rightly the most famous and better studied 
and provides the cleanest example of halo nuclei to date ||5|- [ri|] . 

In halo nuclei, some of the constituent neutrons or protons venture beyond the drop's 
surface and form a misty cloud or halo. Not surprisingly, these extended nuclei behave very 
differently from ordinary ("normal") nuclei lying along the so-called stability valley in the 
chart of nuclides. In particular, they are larger than normal nuclei of the same mass number, 
and they interact with them with larger cross sections as well. In the case of ^^Li, the last 
two neutrons are very weakly bound. Consequently, these neutrons need very little energy 
to move away from the nucleus. There they can remain in their "stratospheric" orbits, 
spreading out and forming a tenuous halo. If one neutron is taken away from ^^Li, a second 
neutron will come out immediately, leaving behind the core of the system, the ordinary 
nucleus ^Li. This result testifies to the fact that pairing, the attraction correlating pairs 
between the least bound particles in a system, plays a central role in the stability of ^^Li. 

It is well known that pairing can radically affect the properties of a many-body system. 
In metals, pairing between the electrons gives rise to superconductivity |T^. In a dilute 



neutron gas, pairing can influence the properties of neutron stars, cold remnants of fierce 
supernova explosions |[T^. Acting on a liquid made out of atoms of the lighter isotope of 
helium (2Hei) it leads to superfluidity [|I4|. In nuclei, it controls almost every aspect of 
nuclear structure close to the ground state and determines, to a large extent, which nuclei 
are stable and which are not @|TB . 



The basic experimental facts which characterize ^^Li and which are also of particular 



relevance in connection with pairing in this system are: a) gLig and s^Lig are stable, s^Liy 
is not, b) the two-neutron separation energy in ^^Li is only 5*2^=0. 294±0. 03 MeV 0] as 
compared with values of 10 to 30 MeV in stable nuclei, c) ^°Li displays s- and p-wave 
resonances at low energy, their centroids lying within the energy range 0.1-0.25 MeV and 0.5- 
0.6 MeV respectively |T§, d) the mean square radius of ^^Li, {r'^Yl'^=?>.hh± 0.10 fm 
is very large as compared to the value 2.32±0.02 fm of the ^Li core, and testifies to the fact 
that the neutron halo must have a large radius (~6-7 fm), e) the momentum distribution of 
the halo neutrons is found to be exceedingly narrow, its FWHM being equal to cr_L = 48 ± 10 
MeV/c for the (perpendicular) distribution observed in the case of the break up of ^^Li 
on ^^C, a value which is of the order of one fifth of that measured during the break up of 
normal nuclei f) the ground state of ^^Li is a mixture of configurations where the two 

halo nucleons move around the ^Li core in s^— and p^— configurations with almost equal 



weight [20,21 



Before discussing the sources of pairing correlations in ^^Li, we shall study the single- 
particle resonant spectrum of ^°Li. The basis of (bare) single-particle states used was de- 
termined by calculating the eigenvalues and eigenfunctions of a nucleon moving in a Saxon- 
Woods potential with spin-orbit and symmetry term (cf. 0, Vol. I, Eqs. (2.281) and 
(2.282)). The continuum states of this potential were calculated by solving the problem in 
a box of radius equal to 40 fm, chosen so as to make the results associated with ^°Li and 
^^Li discussed below, stable. While mean field theory predicts the orbital to be lower 
than the Si/2 orbital (cf. Fig. 1, 1(a)), experimentally the situation is reversed. Similar 
parity inversions have been observed in other isotones of 3°Li7, like e.g. 4^Be7. Shell model 
calculations testify to the fact that the effect of core excitation, in particular of quadrupole 
type, play a central role in this inversion [^,(cf. also ||2^). In keeping with this result. 



we have studied the effect the coupling of the pi/2 and Si/2 orbitals of ^°Li to monopole-, 
dipole- and quadrupole- vibrations of the ^Li core has on the properties of the 1/2"*" and l/2~ 
states of this system. In this study we have used vibrational states calculated by diagonaliz- 
ing, in the random phase approximation (RPA), a multipole-multipole separable interaction 



taking into account the contributions arising from the excitation of particles into the con- 
tinuum states. The couphng strengths were adjusted so as to reproduce the position and 
transition probabihties of vibrational states in the normal nucleus ^°Be |2^. Unperturbed 
particle-hole excitations up to 70 MeV have been included in the calculations and phonon 
states up to 50 MeV have been considered. Within this space there are of the order of 10^ 
dipole states and about the same amount of quadrupole states, exhausting the associated 
energy-weighted sum rules. The resulting vibrational states were coupled to the particle 
states making use of the associated transition densities (formfactors) and particle- vibration 
coupling strengths. A Skyrme-type effective interaction (SLy4) was used to calculate the 
monopole linear response. The corresponding solutions were obtained in coordinate space 
making use of a mesh extending up to a radius of 80 fm. The monopole response exhausts 
94% of the EWSR considering the summed contributions up to 40 MeV of excitation en- 
ergy. This response function was discretized in bins of 300 keV and the resulting states 
coupled to the single-particle states making use of the corresponding transition densities 
and particle-vibration coupling strengths. Similar calculations were performed to determine 
the properties of the L=0,1 and 2-vibrational states of ^^Li needed in connection with the 
discussion of the pairing induced interaction carried out below. In this case the strength 
of the separable dipole interactions has been slightly changed from the value used in the 
previous calculation so as to provide an overall account of the experimental dipole response 
in ^^Li p5|j26[| . For the quadrupole strength the same value used before was employed. Be- 
cause the calculations have been carried out on the physical (correlated ) ^^Li ground state, 
the transitions associated with the vibrational states involving the pi/2 and the Si/2 states 
have been renormalized with the corresponding occupation numbers resulting from the full 
diagonalization. 

In the study of the single-particle resonances of ^°Li we have considered not only the 
particle- vibration coupling vertices associated with effective-mass-like diagrams (upper-part 
graph of Fig. 1, 1(b)) and leading to attractive (negative) contributions to the single- 
particle energies, but also those couplings leading to Pauli principle (repulsive) correction 



processes associated with diagrams containing two-particles, one-hole and a vibration in the 
intermediate states (lower-part diagram of Fig. 1, 1(b)). Because of such Pauli correction 
processes, the pi/2 state experiences an upward shift in energy, arising from the coupling 
of this orbital to the hole-state through quadrupole vibrational states, in keeping with 
the fact that the (P1/2P3/2) particle-hole excitation constitutes an important component of 
the quadrupole vibration wavefunction. As a consequence, the pi/2 state becomes unbound, 
turning into a low- lying resonance with centroid Eres ~ 0.5 MeV. Due to the coupling to 
the vibrations the s— states are instead shifted downwards. In fact, in this case there are 
essentially no (repulsive) contributions arising from the Pauli correction processes. On the 
other hand, (attractive) effective-mass-like processes with intermediate states consisting of 
one particle plus a vibrational state of the type (si/2 x 0^), {pi/2 x 1 ) and x 2+) lead to 
a virtual state with E^i^t = 0.2 MeV (cf. Fig. 1, 1(b)), the coupling to quadrupole vibrations 
providing by far the largest contribution to the total energy shift. The important difference 
between the distribution of the single-particle strength associated with the resonant state 
Pi/2 and the virtual state Si/2 can be observed in Fig. 1, 1(c), where the partial cross section 
ai for neutron elastic scattering off ^Li is shown. While cTp displays a clear peak at 0.5 MeV, 
as is a smoothly decreasing function of the energy. On the other hand, a small increase in 
the depth of the potential felt by the s— neutron will lead to a (slightly) bound state, hence 
the name of virtual |^ . 



At the basis of the variety of pairing phenomena observed in systems so apparently 
different as atomic nuclei and metals is the formation of Cooper pairs p8| , p9| . In the case of 



metals. Cooper pairs arise from the combined effect of Pauli principle and of the exchange 
of lattice vibrations (phonons) between pairs of electrons (fermions) moving in time reversal 



states lying close to the Fermi energy. In the superconducting or BCS state of metals [jT2 
Cooper pairs are strongly overlapping. In fact, there are on average 10^ pairs which have 
their centers of mass falling within the extent of a given pair wavefunction. In spite of the 
modest number of Cooper pairs present in the ground state of atomic nuclei (< 10), BCS 
theory gives a quite accurate description of pairing in nuclei the analogy between 



the pairing gap typical of metallic superconductors and of atomic nuclei being very much 
to the point ||31]]. On the other hand, the finiteness of the nucleus introduces in the BCS 
treatment of pairing important modifications (quantal size effects (QSE), ||30|, [0- [0). 
In particular, while in the infinite system the existence of a bound state of the (Cooper) 
pair happens for an arbitrarily weak interaction PSI, in the nuclear case this phenomenon 



takes place only if the strength of the nucleon-nucleon potential is larger than a critical 
value connected with the discreteness of the nuclear spectrum. In fact, calculations carried 
out making use of a particularly successful parametrization of the (bare) potential (Argonne 
potential ||35|), testify to the fact that the nuclear forces are able to bind Cooper pairs in 
open shell nuclei (leading to sizable pairing gaps (1-2 MeV) |3^), but not in closed shell 
nuclei, the most important contributions to the nucleon-nucleon (pairing) interaction arising 



from high multipole components of the force |30 



The situation is however quite different for the "open shell" nucleus ^^Li. In fact, allow- 
ing the two neutrons to interact through the Argonne potential produces almost no mixing 
between s— and p— waves, but essentially it only shifts the energy of the unperturbed (res- 
onant) states 5^^2(0) ^1/2(0) tiy about 80 keV without giving rise to a bound system. 
Similarly, making use of the matrix elements of the same nucleon-nucleon potential in con- 
nection with the BCS equations does not lead to a solution but the trivial one of zero pairing 
gap (A^ = 0, U^Vi, = 0) At the basis of this negative result is the fact that the most 
important single-particle states allowed to the halo neutrons of ^^Li to correlate are the 
Si/2, Pi/2 and orbitals. Consequently the two neutrons are not able, in this low-angular 
momentum phase space, to profit fully from the strong force-pairing interaction, as only the 
components of multipolarity L = 0, 1 and 2 of this force are effective in ^^Li because of 
angular momentum and parity conservation . 



In keeping with this result and with those of ref. ||4^, and making use of the fact that 
^^Li displays low-lying collective vibrations p5|,|26|, one can posit that the exchange of these 



vibrations between the two neutrons is the main source of pairing available to them to 
correlate. In fact, allowing the two neutrons to both exchange phonons (induced interaction. 



Fig. 1, 11(a)), as well as to emit and later reabsorb them ( self-energy correction, Fig. 1, 1(b)), 
leads to a bound (Cooper) pair, the lowest eigenstate of the associated secular matrix being 
Egs = -0.270 MeV. Adding to the induced interaction the nucleon-nucleon Argonne potential 
one obtains Egg = -0.330 MeV, and thus a two-neutron separation energy quite close to the 
experimental value. Measured from the unperturbed energy of a pair of neutrons in the 
lowest state calculated for ^°Li, namely the s-resonance {Eunp = '^Eg^^^ = 400 keV, cf. Fig. 
l,I(b)), it leads to a pairing correlation energy Eo = Eunp — Eg ^, = 0.730 MeV (cf. Fig. 1, 
11(b)). From the associated two-particle ground state wavefunction \E'o(n, '^2)(= (^i, "^210^)), 
one obtains a momentum distribution (whose FWHM is a± = 56 MeV/c, for ^^Li on ^^C) and 
ground state occupation probabilities of the two-particle states 3^^2(0)) ^1/2(0) '^5/2(0) 
(0.40, 0.58 and 0.02 respectively, cf. Fig. 1, 11(b)) which provide an overall account of the 
experimental findings. The radius of the associated single-particle distribution is 7.1 fm. 
Adding to this density that of the core nucleons one obtains the total density of ^^Li. The 
associated mean square radius (3.9 fm) is slightly larger than the experimental value. 

The spatial structure of the Cooper pair described by the wavefunction \l/o(^i) ''^2) is 
displayed in Fig. 2. The mean square radius of the center of mass of the two neutrons is 
(r^^)^/^ = 5.4 fm. This result testifies to the importance the correlations have in collecting 
the small (enhanced) amplitudes of the uncorrelated two-particle configuration 5^^2(0) 
the region between 4 to 5 fm, region in which the ^^^2(0)? helped by the centrifugal barrier, 
displays a somewhat larger concentration (Fig. 3). From the above results, it emerges that 
the exchange of vibrations between the least bound neutrons leads to a (density-dependent) 
pairing interaction acting essentially only outside the core (cf. also ref. [0]). To be noted 
that the long wavelength behaviour of these vibrations is connected with the excitation 
of the neutron halo, the large size of which not only makes the system easily polarizable 
but provides also the elastic medium through which the loosely bound neutrons exchange 



vibrations with each other M 



The average mean square distance between the halo neutrons is < r^2 >^^^~ 9.2 fm, in 
keeping with the fact that the coherence length [|12| associated with Cooper pairs in nuclei is 



larger than the nuclear dimensions thus preventing the possibility of a nuclear supercurrent 
(cf. e.g. [0 Vol. 11, p. 398). On the other hand, this value of < > does not prevent 
the two correlated neutrons to be close together. The corresponding (small) probability (cf. 
Fig. 2) being much larger than that associated with the uncorrelated neutrons (cf. Fig. 3). 

Similar results as those reported above are obtained solving the BCS equation for the 
two-neutron system making use of the matrix elements used in the diagonalization, sum of 
those of the nucleon-nucleon Argonne potential and those of the induced interaction. In 
this case, the correlation energy is Eo= 0.7 MeV, the separation energy of the two neu- 
trons becoming S2n = 0.3 MeV. The radial structure of the projected BCS wavefunctions 
J2u>o(yu/Uu)fu{^i)ip^{f2) displays a spatial structure quite similar to \t'o(ri,r2), the admix- 
ture of s- , p- and d- two particle configurations being now 46% and 51% and 3% respectively. 
The coherence length ^, that is the mean square distance between the two neutrons forming 
the Cooper pair, is in this case, (^"12)^''^ = 7.8 fm. 

Arguably, the understanding of halo nuclei is the single, most important issue of nuclear 
structure research still awaiting a satisfactory explanation. We have shown that a substantial 
advance towards this goal is made by properly characterizing the role the surface of the 
system plays in renormalizing the bare nucleon-nucleon potential and the single-particle 
motion of the nucleons. In fact, we find compelling evidence which testifies to the fact that 
the mechanism which is at the basis of the presence of a low density halo in ^^Li, is the 
coupling of weakly bound nucleons to long wavelength vibrations of the system leading to 
pairing instability and thus to a bound system. This result suggests a general strategy for 
designing nuclei with a large excess of neutrons or protons which may prove valuable in the 
current exploration of the drip line of the chart of nuclides, and thus in the design of new 
nucleon species: select those systems which display, under an increase of the excess of one 
type of nucleons, a marked softening of the long wavelength linear response. It is likely 
that such systems could venture far inside the region delimited by the drip lines, in keeping 
with the fact that the associated vibrational modes are expected to be an important source 
of (induced) pairing interaction, and thus to contribute significantly to the stability of the 
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system. 
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Caption to the figures 
Fig. 1 

(I) Single-particle neutron resonances in ^°Li. In (a) the position of the levels 
Si/2 and pi/2 calculated making use of mean field theory is shown (hatched area and thin 
horizontal line respectively). The coupling of a single- neutron (upward pointing arrowed 
line) to a vibration (wavy line) calculated making use of the Feynman diagrams displayed 
in (b) (schematically depicted also in terms of either solid dots (neutron) or open circles 
(neutron hole) moving in a single-particle level around or in the ^Li core (hatched area)), 
leads to conspicuous shifts in the energy centroid of the Si/2 and pi/2 resonances (shown by 
thick horizontal lines) and eventually to an inversion in their sequence. In (c) we show the 
calculated partial cross section ui for neutron elastic scattering off ^Li. 

(II) The two-neutron system ^^Li. We show in (a) the mean-field picture of ^^Li, 
where two neutrons (solid dots) move in time-reversal states around the core ^Li (hatched 
area) in the si/2 resonance leading to an unbound 5^^3(0) state where the two neutrons are 
coupled to zero angular momentum. The exchange of vibrations between the two neutrons 
shown in the upper part of the figure leads to a density dependent interaction which, added to 
the nucleon-nucleon interaction, correlates the two-neutron system leading to a bound state 
|0"^), where the two neutrons move with probabihty 0.40, 0.58 and 0.02 in the two-particle 
configurations 5^^2(0)) ^1/2(0) (^5/2(0) respectively. 

Fig. 2 

Spatial structure of two-neutron Cooper pair. The modulus squared wavefunction 
|*o(rL,^2)p = |(ri,^2|0+)p (cf. Fig. 1, II (b)) describing the motion of the two halo 
neutrons around the ^Li core (normalized to unity and multiplied by IGn^rfr'^) is displayed 
as a function of the cartesian coordinates X2 — r2COs(^i2) and 1/2 — r2sin(^i2) of particle 2, 
for fixed value of the position of particle 1 {ri—2.5, 5, 7.5 fm) represented in the right panels 
by a solid dot, while the core '-"Li is shown as a red circle. The numbers appearing on the 
z-axis of the three-dimensional plots displayed on the left side of the figure are in units of 
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Fig. 3 Spatial distribution of the pure two-peirticle configurations s\i^ (0) 

and ^^^2(0) as a function of the x- and y-coordinates of particle 2, for a fixed value of the 
coordinate of particle 1 (ri=5 fm). For more details cf. caption to Fig. 2. 
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